Topology of event horizon in axially symmetric spacetime: 
Classification of Maxwell set with symmetry 



The final state of the black hole is well understood and is simply described by at most three parameters. 
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Abstract 

The crease set of an event horizon is studied in a spacetime with discrete or continuous 
symmetry. It determines possible topologies on spatial sections of an event horizon. 
We thereby investigate the classification of stable topological structure of the crease 
sets in a spacetime with any symmetry. In practice, we show the classification of the 
crease set in axially symmetric spacetime. By that we realize the topological structure 
of axially symmetric event horizons. We will fine that many new topological structures 
^ ' become stable which is not stable without symmetry. 

§ ! 1 Introduction 
in 

o 
in 

On the other hand, an early stage of black hole formation is not well understood since there is much 
O . variety of their appearance. 

In the topological viewpoint, the final state was investigated by many authors and now it is known 

5h . 

bJ) 1 that under some reasonable conditions such as asymptotic flatness and the weak energy condition, each 

component of the black hole region is topologically trivial, i.e. simply connectedpQ. On the other hand, 
there were numerical simulations which suggest non-trivial topologies of event horizons in an early stage 
of black hole formation 2; .3- . There has been some confusion, but it is now well understood that even 
though the black hole region in the spacetime is simply connected, there are many possible topologies of 
spatial sections. 

In particular, it is revealed that various topologies of the spatial sections of black hole are determined 
by the endpoint set, or similarly, the crease set, of the event horizon by the author^]. Therefore the 
crease set of an event horizon, is an important object which is independent from the choice of time slices 
and which determines qualitative properties of the event horizon. Thus, in order to restrict the physically 
possible topologies of the black hole, it is important to restrict the possible (stable) structure of the 
crease set. Since the crease set can be regarded as a singularity of a distant function (i.e. the Maxwell 
set) determining the generators of the event horizon. The singularity theory of real maps can classify the 
stable topological structure of the crease set locally |Sj 00 • 

These possible topological structures should be studied in numerical experiments of gravitational 
collapse 2 '■'> s]. Nevertheless, it might be very hard since it requires a huge data set of a geometry in the 
entirety of numerically generated spacetime, to determine its event horizon. To reduce the cost of the 
numerical calculation, it may be successful to impose any symmetry on a black hole spacetime. Indeed it 
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is in axially symmetric gravitational collapse that some remarkable results are reported about non-trivial 
topologies of the event horizon[2][3]. 

In the symmetric spacetime, however, is the previous resultjS] how to be applied to? In a rigorously 
symmetric spacetime, their perturbation should be also symmetric. Then the stability of the function 
determining the generators of an event horizon will be altered. We need to study the classifications of the 
crease set again in this symmetric situation. By this, some new crease sets are added to the classification 
without symmetry. 

This investigation will become relevant to also an almost symmetric realistic spacetime, if the sym- 
metry is dynamically stable and if non-symmetric perturbation is always sufficiently small. In such a 
situation, the crease set in the almost symmetric spacetime will be approximated by the symmetric clas- 
sification of the crease set in rigorously symmetric spacetime, by neglecting very fine structures of the 
crease set (c.f. [§]). 

In the present article, we investigate the classification of the Maxwell set under discrete and continuous 
(especially, reflection and axial) symmetry. Of course, the classification includes the event horizons 
reported in references 00. 

In the next section, we mention the definition of variational principle to determine the null generators 
of event horizons in symmetric spacetime. The crease set is mathematically related to the Maxwell set 
of a potential function in the third section. We investigate the classification of the symmetric Maxwell 
set in the fourth section. The fifth section gives a concrete classification in reflection symmetric three- 
dimensional spacetime and axially symmetric four-dimensional spacetime. Then we will discuss the spatial 
topology of axially symmetric black hole. The final section is devoted to summary and discussions. 

2 Fermat potential in symmetric spacetime 

An event horizon is generated by null geodesies. A future event horizon TL cannot have future endpoints, 
but can have past endpoints if it is not eternal. As is pointed out in 0], the endpoint set £ of a horizon is 
an arc-wise connected acausal set. Points u G £ are classified by the multiplicity m(u) of u, the number 
of the generators emanating from u: 

C := {u e £ | m(u) > 1}, 

V := {u e £ | m{u) = 1}. (1) 

The set C is called the crease set of the horizon. The crease set contains the interior of the endpoint set, 
i.e., the closure of C contains £ |11|. The crease set C equals the set of points of £ on which the horizon 
is not differentiable, i.e., the horizon is diffcrentiable at u € £ if and only if u € T> [111 I12j . 

The horizon TL is the envelope of the light cone starting from the crease set C which is an arc-wise 
connected acausal subset of TL. If the spatial section of the horizon is a topological sphere at late times, 
the topology of the spatial section of the horizon can be nontrivial only at the crease set and the topology 
is completely determined by the time slicing of the crease set, which is studied in Ref. in detail. In 
particular, when the crease set is a single point, each possible spatial section of the horizon topologically 
is a single sphere. On the contrary, when the crease set has a disk-like structure, the horizon can have 
toroidal or higher-genus spatial sections. One would see the coalesce of horizons if the crease set has a 
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line- like structure Therefore by classifying the structure of the crease set, we will know all possible 
topologies of the horizons. Here we do not assume that the spatial section of the horizon in the future is 
a sphere. 

The crease set can be determined by Fermat's principle in simple stationary spacetime[5]. In non- 
stationary spacetime, we can extend Fermat's principle and find a variational principle about light paths, 
imposing some appropriate causality condition such as global hyperbolicity. Here we show an example of 
the construction of the Fermat potential, and see how the symmetry of spacetime impose a condition to 
the Fermat potential. 

Let us assume that the spacetime M. is smooth and is globally hyperbolic from a smooth Cauchy 
surface S which is diffeomorphic to K n . Furthermore, we consider spacetime of gravitational collapse, 
namely, we assume that the event horizon TL is in the future of S. 

By global hyperbolicity, there are always an appropriate smooth global time coordinate t : M — > K 
and a timelike vector field T such that dt(T) = 1. The spacetime M. is foliated by Cauchy surfaces 
St = {q G A4\t(q) = t}. The vector field T — d/dt defines a smooth projection tt from M. into the 

S — S+ n i 



dt 



= {7(*)lSr = T Mto) = q,t ei}. (2) 



Conversely, there is a diffeomorphism 

(f> :Rx S ^ M, 

*((£(*, u)) =t, n(<p(t,u)) = u. (3) 
Because TL is achronal, the restriction of tt on TL is injective and has an inverse, which we denote by ip: 

xp : tt(TL) -» TL, 

ip(u) 6 TL, n(ip(u)) = u. (4) 

The map ip is Lipschitz|10|. 

In the present article, we suppose that the spacetime possesses any symmetry. This means that there 
exists a symmetric Cauchy surfaces Sst and there are any isometries ip acting on each spatial hypersurface 
Sst isometrically, 

<p:(M,g)~(M,g) (5) 
s.t. tp(S S t) = S st (6) 
<p*g = 9,^*1* = h , (7) 

where h is a spatially induced metric on Sst- We consider both cases of discrete symmetry and continuous 
symmetry. For these cases their isometries form discrete group and continuous group, respectively. 

Then we adopt a time function of this symmetric timeslicing as a symmetric global time coordinate 
t(q) and its normal time vector T — d/dt _L Sst defines symmetric smooth projection n. By the fact that 
d/dt and the Killing vector generating the isometry JZJ are orthogonal, t can be a synchronous time, 

d d\ 

dl'di =COnSt (8) 
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We take some (sufficiently large) t = ti and assume that Sgt ± H TL is diffeomorphic to a compact 
manifold M . We consider M as a fixed submanifold embedded in Ss so that iSs tl DTL = tp(M). Consider 
a neighbourhood [/ of ir(TL n J~(Sgt 1 )) in <Ss- For x G M and u £ U we define Fermat potential as 
follows: 

F(x, it) := - sup{f G R\<f>(t, u) G J" (9) 
Obviously if there is an event horizon, it should be invariant under the isometry tp 

tp(H) = TL. (10) 
Then M is also invariant under the isometry tp 

<p(M) = M. (11) 

Furthermore we take such U that it is also invariant subset, tp(U) = U. 

The minimum points x of F corresponds to the generator of TL through it. Our definition Q is the 
generalization of the geodesic distance function to the non-static spacetime. 

From the above construction of the Fermat potential, it will be understood how the Fermat poten- 
tial reflects the symmetry of spacetime. Since we adopt a spatially symmetric timeslicing Ss and its 
synchronous normal time vector d/dt, F(x,u) is invariant under isometric transformation. In this con- 
struction the isometry acts on both state X and control U spaces. Then pull back of the Fermat potential 
satisfies 

tp*F(x, u) = F(<p(x)Mu)) = F ( x > «). (12) 
where the isometry <p acts on M and U as a invariant subset of M. from (|10fl . (|ll|) and 

7T O tf — if O 7T. (13) 

From ©, the crease set C is given by 

C = VK^Maxwcll^)), (14) 

where BMaxwQii(F) is the Maxwell set of F where F has two or more minimum points. A precise mathe- 
matical framework to study its properties is given in [Sj. By (I12II . C is also invariant tp(C) = C. 

The goal of our present study is to classify all generally possible structure about the singularities 
(Maxwell set) of the Fermat potential determining the event horizon. As shown in ref.jS], the generic 
structure will be given by studying singularities of stable Fermat potentials concretely. Then our main 
task is to give a classification of the stable functions satisfying When one usually discusses the 

bifurcation structure, so-called caustics, of a system we analyze the Fermat potential locally in the context 
of a function germ. However, our main object here is not the bifurcation set but the Maxwell set (the 
difference is clarified in ref.[l)|). The problem is not purely local but is rather semi-local. The definition of 
Maxwell set is local in control (parameter) space U but is non-local in state (variable) space M. To treat 
this, one introduces function multigerms and classifies stable multiunfoldings and their Maxwell setsjS]. 
On the other hand, our Fermat potential in the present study is restricted to a symmetric one ()12ll . To 
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give a symmetric multiunfolding, non-local [/-variable is prepared as an orbit of the isometry and also 
non-local Af-variable to include non-isolated critical points of a symmetric Fermat potential. In the next 
section, we give these concepts as a G-orbital multigerm of unfolding. 

3 mathematical preparation 

The mathematical framework to study the Maxwell set was given in [S] . There the semi-local characteristic 
of a Fermat potential is investigated in the context of a multigerm. Nevertheless it needs some additional 
preparation to study the Maxwell set with any symmetry. It is caused mainly by the change of a function 
space of perturbation. That requires the extension of a multigerm so as to represent an isometry orbit 
of a multigerm and competence of critical points mapped by the isometries. 

First of all, we prepare the same notions as in JjjJ. A function unfolding F : M x U — > K can be 
considered as a family of functions f u : M — ► K with /„ = F(», u). The Maxwell set of a function unfolding 
is the set of all values of the parameters for which the minimum is attained either at a non-Morse critical 
point or at two or more critical points. 

Definition 1 (Maxwell set). For a function unfolding F : M x U — > K on a compact manifold M , the 
Maxwell set -BMaxweii(-F) of F is a subset ofU given by 

Byiaxwc\\{F) :— {u S U\f u has two or more global minimum points}. (15) 

In the following we sometimes make x as a representative of the state variables and u the control 
variables. Of course, F is not a global function on 1" x R r rather a function on manifold M x U, where 
M is a compact manifold and U is an open subset of K r . 

In the investigation of the Maxwell set we mainly focus on its local structure because the global 
structure is obtained by the combinations of local ones. Below we extensively use the notion of the germs 
of objects which provides the best way to characterize their local structure. The rigorous definition of 
germ is well known and will be given elsewhere[H|. Let M, N be C°°-manifolds. We denote the set of 
C°°-maps from M to N by C°°(M,N). 

The definition of the Maxwell set requires the global information of the function unfolding F. A 
simple but crucial observation is, however, that to determine the local structure of the Maxwell set, 
i.e., the Maxwell set germs, we only need the local information of F around its global minimum points 
Px, ...,Pk 6 M. We generalize the notion of germs to that of multigerms. Let be a fc-tuple of distinct 
points of M , i.e., 

:={( Xl ,...,x k ) &M k \ Xi ^ Xj for Mi}. (16) 

Definition 2 (Multigerm). Let (x u ...,x k ) e A'F k \ A fc-fold map germ / : (M, (x\, ...,Xk)) — > N, or 

[f]xs. x k , is the equivalence class of f G C°°{M,N), where two maps are equivalent if they coincide on 

some open subset of M which contains x\, ...,Xk- A fc-fold unfolding germ F : [M x U, ((xi, Xk), Uo)) — » 
N, or [-P'](a: 1 ,...,x ii ),Mo' * s ^ e equivalence class of F G C°°(M,M.) where two functions are equivalent if they 
coincide on some open subset of M x U which contains (xx, uq), (xk, uo). A k-fold germ is also called 
as a multigerm. 
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A multigerm can be considered as A:-tuple of simple germs. For example, a map multigerm [-P'](xi,...,x fe ),«c 
can be considered as fc-tuple of function germs ([/i]xu •••) [fk]x k )- 

To study the Maxwell set under any symmetry, first we prepare a set of multigerms they are generated 
by a transformation group. 

Definition 3 (G-orbital multigerm of unfolding). Consider a group of transformation G. Let 
(x u ...,x k ) £ M (fc) . G(((.T 1 ,...,x fc ),u )) - {9{{{x u ...,x k ),u ))\g £ G} C AfW x U is the orbit of G 
through ((xi, x k ), uq). A G-orbital multigerm of unfolding Fq '■ (M x U, (G((xi, x k ), Uo))) — > N, or 
[Fg] G{(xi,...,x k ),u )> * s the equivalence class of F £ C°°{M x f7,R.) where two functions are equivalent if 
they coincide on some open subset of x [7 which contains G(((xi, x k ),Uo)). 

A G-orbital multigerm can be considered as a suite of multigerms. For example, if G is a discrete 
group a G-orbital multigerm [FG]G((x 1 ,...,x k ),u ) can be considered as multituple of k-tuple function germs 
(([flU,..., [fl] Xk ), (iflU Xl ), [#],(«»)). (•••), ■■■) where g £ G. 

If a representative of Fq is invariant under actions of g £ G, by definition its G-orbital multigerm 
of unfolding is invariant under actions of G. Here it should be noticed that discrete G gives point-wise 
equivalence similar to usual multigerm otherwise the equivalence is not point-wise. It is also commented 
that when Xi or uo is a fixed point of G the equivalence becomes point- wise there even if G is continuous. 

It is required to define not point-wise equivalence of maps to discuss the diffeomorphism class of this 
G-orbital multigerm of unfolding. 

Definition 4 (Map germ at subset). Subset S C M is given. Maps /, g £ C°° (M, N) are equivalent 
at subset S if there is a neighbourhood W of S such that f\w = g\w- A map germ / at S, [f]s, is the 
equivalence class of f . It is also denoted by f : (M, S) — > N or f : (A/, S) — > {N, f(S)). 

Henceforth we call germ at subset simply germ and sometimes also germ is omitted. Examples of 
map germs at subset include function germs and diffeomorphism germs. 

As easily seen, that G-orbital multigerm of unfolding is not local also in U. When we find the Maxwell 
set of it, it will be expressed as a non-local equivalence class of set. 

Definition 5 (Set germ at subset). Subset S C U is given. Subsets X, Y ofU are equivalent at subset 
S ifXdW — Y(~\W for each neighbourhood W of S . A set germ {X, S) of X at S is the equivalence 
class of X. Set germs (X, S) and (Y, S') are diffeomorphic, (X, S) ~ (Y, S'), if there is a diffeomorphism 
germ <j) : (M, S) -> M such that (<j)(X), <j>(S)) = (Y, S'). 

Next we give a fundamental concepts to investigate the stability of unfolding under symmetry |14j . 

Definition 6 (Symmetry preserving diffeomorphism). Let (A4,g) be a symmetric space with isom- 
etry group G. A diffeomorphism <f> on (M,g) is a symmetry preserving diffeomorphism (SPD) if it keep 
the action a £ G on (A4,g) isometry. If a diffeomorphism (ft is an SPD then (ft : a ^— > (ft o a o <fr x is an 
automorphism of G. In the following, G mainly acts on invariant subset M x U as 

tp: M xU h-> M xU (17) 
(x,u) £ M x U,<p(x,u) = (ip(x),ip(u)). (18) 

Then SPD is sometimes considered as the diffeomorphism on M x U besides M and U satisfying above 
condition. 
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In the present article we discuss the stability defined by this SPD, since the subset of the function 
space we treat is invariant under the SPD. 

Definition 7 (Right SPD equivalence). Function germs f : (M, x) — > R and g : (M, y) — > R are right 
equivalent, [/] x ~ [g] y , if there exists a SPD germ (f> : (M, x) — > (M, y) and a € R smc/i i/iat f = go(j) + a 
holds as an equality of function germs at x. 

Definition 8 (Right SPD equivalence at the minimum points). Unfoldings F : (M x U,M x 
{u }) — > R and G : (M x U,M x {v }) — ► R are right equivalent at the minimum points 
i/ ifte following conditions hold: 

(1) The functions f Uo — F(»,u n ) and g Vg = G(», vo) have the same number of global minimum points, 
pi,...,Pk and qi, ...,qk, respectively. 

(2) There exists a SPD multigerm 4> : (M x U,((pi, ...,pk),Uo)) — ► (M x U, ((qi, ...,Pk), vo)), a SPD 
germ ip : (U, uq) — > (U, vq), and a function germ a : (U, uq) — > R such that 

F{x,u)=G{4>{x,u),i){u))+a{u) (19) 

holds with both sides being function multigerms at ((pi, ...,Pfc), Uo). 

The Maxwell set germ of an unfolding is determined only by the unfolding multigerm at the minimum 
points: 

Proposition 1. If unfoldings F : (MxU,Mx {u Q }) — > R and G : (M xU,Mx {v Q }) -> R are right SPD 
equivalent at the minimum points, then their Maxwell set germs are equivalent under SPD transformation. 



Proof. Follows directly from the definitions of right SPD equivalence and of Maxwell set germs. □ 

Below, we will determine the topological structure of C°° (M, N) where all the maps that we treat are 
included. We define a topology of C°°(M, N) by the r-jet space J r (M, N) below. 

Definition 9 (Jet space). Let f G C°°{M,N). The r-jet j r f{xo) of f at x is the equivalence class of 
f in C°°(M,N) where two maps are equivalent if all of their s-th partial derivatives with 1 < s < r, in 
some coordinate systems of M and N, coincide. The r-jet space of C co (M,N) is defined by 

J r {M,N) := {ff(x Q )\f e C°°(M,N)}. (20) 

(tn ~\~ t \ 
I -dimensional manifold, where m = dimM and 

n = dim TV. 

Now we endow the space C°°(M,N) with the Whitney C°° topology. 
Definition 10 (Whitney C°° topology). For an open subset O of J r (M,N), let 

W r (0) := {/ £ C°°(M,N)\ff(M) C O}. (21) 
The Whitney C°°topology on C oc (M,N) is the topology whose basis is 

oo 

\J{W r (0)\0 is an open subset ofJ r (M,N)}. (22) 

r=0 
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Hereafter we treat that C°°{M,N) as a topological space with the Whitney G°° topology. Now we 
can define stability of the Maxwell set using this topology. 

The concepts of the jet space and the Whitney G°° topology are used as in the ref.|S]- Especially, 
multiunfolding was investigated by the jet space for C°°(M x U,M). In the present article, however, it is 
not this C°° (M x U, R) in which we investigate the stability of a G-orbital multigerm of unfolding but 
its subset whose elements are invariant under an isometry group G (112(1 . We call it G-C°° (M x U, R). 
From (|12f) . it can be identified with G°°(^p,R). Here we propose that the topology of the subset is 
relative topology and equivalent to the topology of C°° ( M q , R) in the context of the Whitney topology. 
G-G°°(M x 17, R) is invariant subset under the pull back of SPD. 

Now, we have completed the preparation. We will investigate stable structure of the function unfolding 
and its Maxwell set. 

Definition 11 (G-stable Maxwell set germ). A G-orbital multigerm of unfolding, 



is stable with respect to the Maxwell set if for each neighbourhood W of G(uq) there exists a neigh- 
bourhood U of Fq in G-G°°(M x U, R) such that for each F' G G U there exists Vq G W such that 
(-BMaxwellC-Fb)) G(uo)) and (Bwiaxwea(FQ),G(va)) are symmetry preserving diffeomorphic about G. 
We call (BMaxwe\i(F),G(uo)) a G-stable Maxwell set germ. 

Definition 12 (G-stability at the minimum points). A G-orbital multigerm of unfolding Fq : 
(M x U, ({x\, Xk), uo)) — > R is G-stable at the minimum points if for each neighbourhood W of G(uq) 
there exists a neighbourhood U of Fq in G-C°°(M x U, RJ such that for each F' G 6 U there exists 
G(«o) G W such that [-F]g(u ) an( ^ [F']g(v ) are fight SPD equivalent at the minimum points. 

From Proposition ^ we immediately have the following proposition. 

Proposition 2. If a G-orbital multigerm of unfolding Fq ■ (M x U,M x G(uq)) — + R is G-stable at the 
minimum points then (i?Maxweii(-F), G(uo)) is a G-stable Maxwell set germ. 

In this sense, our aim is to classify the G-stable Maxwell set germs of the Fermat potential F. In 
the present article, we will not prove this stability in practice. We simply suggest the stability from the 
result of Ref.[H] without any symmetry. 

To discuss the stability and classification, we study the orbit of diffcomorphism on some standard 
functions in J r (M, R) and its transversality. We will stratify J r (M, R), i.e., decompose the J r (M, R) 
into the union of submanifolds (strata) ^H] ■ 

Definition 13 (Strata). 



F G : (M x U,G((( Xl ,...,x k ),u ))) -> M, 



(23) 



An := {ff(p) G J r {M,R) \ f is regular atp}, 
A k := {ff(p) G J r (M,R) | [f] p ~ [±^ +1 ] }. 



(24) 
(25) 



The following is well discussed |18| : 



Lemma 1. (1) Aq is an open subset of J r (M, R) hence is a submanifold of codimension 0. 

(2) A k is a submanifold of J r (M,M) of codimension k. 

(3) £ := J r (M, M) — (Jfc=o Ah is the union of a finite number of submanifolds of codimension 4 or 
greater: 

£ = WiU...UW S) (26) 

where the codimension is determined for the case M is one- dimensional. It is different from that of 
Ref.jg. 

Definition 14 (Natural stratification). The natural stratification of J r {M, M), where r > 3, is the 
one given by 

S(J r {M,R)) := {A ,A 1 ,...,A 3 ,W 1 ,...,W a }. (27) 

This natural stratification corresponds to the classification of the stable unfolding under the aid of 
well-known transversality theorem. To discuss the Maxwell set we naturally extend this concept to the 
one at several minimum points. For that, we extend the jet space to multijet space. 

Definition 15 (Multijet space). The fc-fold r-jet, or simply, r-multijet, kj r f of f : M — » K at 
(xi, ...,x k ) £ MW is 

kff(xi,...,x k ) -(ffM-'ffW)- (28) 
The r-multijet space feJ r (M, N) is given by 

k J r (M,N) := { k ff(x u ...,x k ) G (J r (M,N)) k \ (x u x k ) G M«}, (29) 
The map k j r f '■ — > k J r (M,N) is called an r-multijet section. 

Definition 16 (Natural stratification of k J r {M, K)). T/ie natural stratification of k J r (M,M) is the 
one given by 

S( k J r (M,R)) :={A fc nXi x ... x X k \X u ...,X k G 5(J r (M,R))}, (30) 

w/iere 

A fc :={(i r /i(Pi).-.i r A(Pfc)) e kJ r (M, R)|/i(pi) = ... = f k ( Pk )}. (31) 

Hence we find that the stable multiunfolding germ is given by multituple of well known single stable 
unfolding under the aid of Mather's multitransversality theorem^!]- This does not change in the case 
with discrete symmetry since the multiunfolding is point-wise. With continuous symmetry we cannot use 
the multitransversality theorem directly since a family of uncountable unfoldings is considered. 

Finally we give a minimum function to investigate the concrete Maxwell set. If the minimum function 
fiF has a singularity at the point uq, then the function f UQ has several global minimal points on the 
manifold M. This will give just the Maxwell set even if the minimum point is not isolated in the case of 
continuous symmetry. 
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Definition 17 (Minimum function). The minimum function fip : U — > R of an unfolding F : MxU — » 
R is given by 

fi F (u) :=rmn{F(x,u)\x € M}. (32) 



There are two cases when the minimum function /ip has a singularity at u: 

(1) the function /„ has several minimum points in M , or 

(2) the number of minimum points changes there. 

In the case (1) u is a point of the Maxwell set. In the case (2) u is not a point of the Maxwell set but 
corresponds to a point of the endpoint set £. 

4 classification of the symmetric Maxwell set 

In the paper[H|, the stable Maxwell set of multhmfolding F : (M x U, M x {uq}) — > R, where M is two 
dimensional manifold and U is three-dimensional open subset, were mathematically classified. There, 
based on well established local investigation of the unfolding 16 , semi-local classification of the universal 
multiunfolding is made and it gives the classification of the Maxwell set. Then using multitransversality 
theorem, it is shown that they are stable. 

On the other hand, some revisions are required for the present study in symmetric spacetime. With 
any symmetry, we should impose a condition for the symmetry to the list of the standard functions 
which are representatives of classification. Then diffcomorphism equivalence of the classification would 
be restricted to SPD equivalence under the request of the symmetry. We realize the restriction by using 
some global coordinates related to the symmetry instead of the local coordinate of neighborhoods. 

That simply results from imposing a symmetry to the Maxwell set but revisions are not only this. 
With any symmetry, we should change the concept of the stability of the unfoldings from general ones 
to symmetric ones fDcf llll Dcf ll2|) since the definitions of perturbation and the stable multiunfolding 
are altered by the SPD and the symmetric subset G-C°°(M x U, R). Consequently, it is expected that a 
new class of multiunfolding will be added to the classification as a stable one. This seriously embarrasses 
the problem. Especially for the case of continuous symmetry, the function space where unfolding and 
perturbation are defined changes its dimensions (strictly speaking, the dimensions of its jet space) by 
infinite dimensions. This implies that an excluded function germ because of its infinite codimensions 
(which means the function germ cannot have universal unfolding) in function space C°°(M x U, R) is 
possible to be with finite codimensions and provide a stable unfolding in G-C°°(M x U,M). To find this 
directly, we should take a systematic survey of all functions in G-C°°(M X U,R) ~ C*°°(^p,R), and 
is not easy because U/G has a boundary at the fixed points of G. In addition to it, because of such 
continuous symmetry, the concept of the usual point-wise germ will fail since the function germ orbits 
along the isometry. Especially at the fixed point of the isometry in U, their minimum points become 
non-isolated. 

In the previous section, a G-orbital multigerm have been defined to resolve the latter problem (Def ■ 
To avoid the former trouble, we first consider a discrete symmetry Gd in low-dimensional spacetime M' , U' 
whose catastrophe is essentially equivalent to the continuous symmetry G c because of U/G c ~ V /Gd- For 
example, the reflection symmetry in (n — l)-dimensional spacetime and axial symmetry in n-dimensional 
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spacetime are so. Since the change from the general function space C°°(M' x f7',R) to the discrete 
symmetric function space Gd-C°°(M' x U',M) is exceedingly small so that one can systematically pick 
up multiunfoldings newly added to the stable classes. 

First, we recall a symmetry condition about the unfolding. We suppose the spacetime (A4,g) is 
symmetric and there exists a group of isometry G 3 ip, ip*g = g (eq.Q). As stated in the section 2, M 
and U are invariant subspace by the construction, 

<p(M) = M, (p(U) = U, ip(M xU) = MxU. (33) 

As the discussion in the section 2, the unfolding is pulled back by tp, and should be invariant under 
the isometric transformation, 



F(x, u) = ip*(F)(x, u) = F(ip(x),(p(u)). (34) 

Therefore, we confirm that the Maxwell set should be invariant by the isometry group. To find a classifi- 
cation of the Maxwell set in this symmetric spacetime, we will list the G-stable multiunfolding. It should 
be a G-orbital multigerm of unfolding (Def 0) . 

4.1 discrete symmetry 

In the case of discrete symmetry investigation is easy to perform. Here, we suppose the isometry group 
is a discrete group Gd 3 tfid- 

4.1.1 Class A multiunfolding 

If the symmetry is discrete, as long as there is no fixed point, the classification of stable multiunfolding 
is essentially same as the case without symmetryj^j. Without fixed points, locally Um * Uu ^ Mm. x u v ^ 
where Um C M and Uu C U is a open subset of M and U, respectively. Then the identification 

kjr ^m^u r ^ „ kJ r (U* x Uu ^ (35) 

implies that the catastrophe of this case is equivalent to the catastrophe without symmetry. 

By the action of the isometry l|34|l on an unfolding F, we generate a G^-invariant G^-orbital multigerm 
of unfolding from a multigerm of unfolding since G^-orbital multigerm is point-wise germ. We write this 
G^-invariant G^-orbital multigerm of unfolding as Gd [F] . 



G d [F] - [Gd[F]]{(( ai! ... iXfc ) j „ ) ! ((^ (i ( Xl ) ; ...^ (i ( !Cfe )) ! ^(„ )) ! ...|^ eG<J } (36) 

= [[ F ]((x u ...,x k ),u ), [ < t ) d F ]{<p d (x 1 ),...,4 )d (x k ))^ d {uo)^ •••] = W: <f>d F > ■•■]' ( 37 ) 
4>d G G d . (38) 

In [Xi, X%, ...], [ ] means that each element Xi does not share uo and is not competent to others, differently 
from (X U X 2 ,...). 

There is no need to care the change of the function space (|35(l . To get a representative function of 
the classification we simply use a function F in the case of the classification without symmetry. 
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Though i*o is mapped to ipd(uo) m (EHJi there is no more degrees of freedom for control parameter 
space U of Gd[-F] than those of the original multiunfolding F. Of course it is obviously G^-stable at the 
minimum points if the original F is stable at the minimum points. The reason is because perturbation is 
considered also in kJ r ( UM Q Uu ,R). We call this a Class A multiunfolding in the following. 

4.1.2 Class B multiunfolding 

Next we consider the case with fixed points of the isometry on M or U in the neighborhood of minimum 
points Gd((pi,P2, ■■•Pk)i uq). Though this point in Af( fe ) or jj [ s the fixed point of only several elements 
gf G G in general, considering the subgroup (gf) = Gf C G generated by them, we suppose the fixed 
point is of all elements of the isometry in the present article. General case would be developed by 
combining the discussion of both cases with and without fixed points. 

If a minimum point pi for uq is a fixed point of tfd, i*o should also be on a fixed point of (fd in our 
discussion. This is required in the case of an event horizon. A couple of uq and a minimum point Pi(u$) 
determines a generator (uo,Pi) of the event horizon. Nevertheless, from the symmetry also (ip(uo) , ip(pi) = 
Pi) is a generator of the same event horizon. This means two generators pass a one point pt on Sstt = 
ip(M). Nevertheless this is not allowed since there is no future cndpoint on the event horizon. Then we 
suppose ito is always a fixed point of the isometry here. 

In this case, we must pay attention to the changes of the function space from C°°(M x U,M) to 
Gd-C°° (M x U,M), where unfolding and perturbation is given since the essential control space U / Gd 
has boundary at the fixed point. To care this problem, it is convenient also to represent a standard 
function by not local coordinates but global coordinates to indicate the symmetry. By the action of the 
isometry to the global coordinate, we observe that the codimension of the stratum for multijet space 
fDef ll5|) is possible to be diminished since uq is common among multituple of multiunfolding for the 
Gd-orbital multigerm of unfolding. When we list the strata for the multijet space, in advance, we should 
take extra strata assuming the decreasing of the control parameter. By this, there is possibility that a 
new multiunfolding becomes stable. 

We find Gd-stable multiunfoldings by the following steps. 

1. One prepares strata for a single-jet space which possess minimum points from Def llHI and give its 
universal unfolding, so that the number of control parameter does not exceed the dimension of U. 
This is well understood and is known to be stable by Thorn's theorem in the context of elementary 
catastrophe. They are written in local coordinates. 

2. List the possible multituples of unfoldings and sum the numbers of their control parameters. By 
the isometric transformation of unfolding (|34J) . some of control parameters are identified in the 
global coordinate for the symmetry. Then we manage the list so that also the sum of the number 
of control parameters does not exceed the dimension of U after such identifications. 

3. Imposing the symmetry relation (|34l) on a possible globalizations of coordinate, we will find a Gd- 
invariant G^-orbital multigerm of unfolding. Of course any of the multituples may not generate 
G^-invariant G^-orbital multigerm of unfolding. 
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4. In the same logic as in jS], it can be shown that the G^-invariant G^-orbital multigerm of unfolding 
constructed above, is the G^-stable at the minimum points, since each unfolding at a minimum 
points is G^-stable unfolding even after imposing the symmetry. This is guaranteed by the fact the 
jet section j r F : M xU i-» J r [M x U,M) transversally crosses the strata also in Gj-G°°(M x U,M.) 
as long as the projection p : C°°(M x U, R) i— > Gd-C°°(M x U, K) maps neither the stratum nor 
the jet section degenerate into a point. It is easily checked that the condition l|34|l never does so in 
most of significant cases. 

Since the G^-orbital multigerm of unfolding found in this way is G^-invariant, it should be provided 
as the one generated by the isometry group Gd- Nevertheless it includes competence at the fixed point 
uq among multiunfoldings generated by each elements ifd G Gd, and is not Gd[F\. Then we write such a 
Gd-invariant G^-orbital multigerm of unfolding as Gd(F); 



F, (j> d F, ... share u and compete to each other. We call this a Class B multiunfolding. 
4.2 continuous symmetry 

Now we write continuous isometry group as ip c G G c . To avoid some difficulties caused by the continuous 
symmetry, here we suppose that the catastrophe of the case of this continuous symmetry is equivalent 
to a catastrophe with any discrete symmetry. Here an essential control parameter space is not U rather 
U/G since the isometry does not allow the unfolding in the direction of the isometry. Then essentially 
the classification of stable multiunfolding with the continuous symmetry is already completed, if such a 
corresponding discrete symmetry really exists and if their catastrophe is equivalent, since in the previous 
subsection, we have investigated the procedure to classify the G^-stable G^-orbital multigerm of unfolding. 

In the next section, it will be shown that there is identification between the control space for a 
reflection symmetry in three-dimensional spacetime and that for an axial symmetry in four-dimensional 
spacetime. Since their catastrophes are equivalent, we only have to translate the classification of the 
multiunfolding with discrete symmetry into that with continuous symmetry. 

4.2.1 Class C multiunfolding 

To find a G c -invariant G c -orbital multigerm of unfolding for the continuous symmetry, it is an easy way 
to generate it by the isometry from the unfolding of discrete symmetry as a orbit of the continuous 
isometry. When there is no fixed point in the neighborhood of ((pi, ...,Pfc), uq), it is possible to generate 
G c -invariant G c -orbital multigerm of unfolding from a G^-invariant G^-orbital multigerm of unfolding by 
straightforward correspondence. 

Here one should be careful to start from a multiunfolding in Class A. To find a G^-invariant G^-orbital 
multigerm of unfolding in Class A, it was generated by the isometry group as Gd[F]. Nevertheless, in 
the continuous symmetry, G c -orbital multigerm suffers overlap of neighborhoods U\ and Hi generated by 
an isometry ip c (Ui) = Ui- Therefore the multigerm F{kp c {x), ip c (u)) may be inconsistent to F(x,u) on 



G d (F) = [Gd(F)]{(( Xu ...,x k ),u ),((Mx 1 ),...,<t, d ^k)),u ),...\<t>4eG d } 

= {[ F ]((x u ...,x k ),u )A4>*d F ](4> d (x 1 ),...,4> d (x k )),u a , ■■■) = (F^dF,...), 
4>d G G d - 



(39) 
(40) 
(41) 
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U\ H142- To resolve this inconsistency it is sufficient that the original unfolding F is invariant along the 
continuous isometry in IA\ . This also is realized by the global coordinate of the continuous symmetry. 

To represent the equivalence class of G^-orbital multigerm, multituple notation was possible in the 
discrete symmetry. On the other hand, in this case of continuous symmetry, we write it as a family of 
uncountable members, 

G C [F] = [G c [F]]^ Va ( Xl ) t _ t9c ^ k ^ tVa ( Uo ^ VceGc y ( 42 ) 
= [{[vl F \{ Va {x 1 ) > ..., Va {x k )) >Vc {u )Wc G G c }] = [{(p$F\(p c £ G c }}. (43) 

Since original F is invariant under the isometry ip*F — F in global coordinate it is sufficient to write as 
[F\Go]. 

Now we need to revive a dummy variable u t of U , which is omitted by the quotient U/G c , to decide 
the Maxwell set in U. Since the generator of isometry d/dt should not effect the catastrophe, that can 
be achieved by adding a Morse function from the splitting lemma Therefore it is sufficient to add 
a St 2 with ip c (8t) = St in order to find G c -invariant F. The dummy variable u t is included in this as 
St = St(u t ). This case is named as a Class C multiunfolding. 

4.2.2 Class D multiunfolding 

At the fixed point, also the G c -invariant G c -orbital multigerm of unfolding is generated by the continuous 
isometry since the fixed point of G c is same as that of Gd- Using global coordinates, it is possible to 
rewrite such multiunfoldings in symmetric form, 

G C (F) = [G c {F)] {{lfa{xi) ^,,^ c ( Xk y huo) \ Vc£Gc] (44) 

= ({[^C^](^a(zi),.MVc(* fe )),«o)kc G G c }) = ({<p* c F\tp c e G c }). (45) 

Since original F is G c -invariant just like the Class C multiunfolding, it is sufficient to write (F ). Here 
(F) should be written as a standard function using a global coordinate in the direction of the continuous 
isometry to clarify that its minimum points make an orbit along the isometry. This will be found by 
discussion in a following concrete situation. Instead of the Morse function, a global function which is 
locally right SPD equivalent to a Morse function are added to revive a dummy variable (see next section) . 

By the splitting lemma[TH| . there is no need to discuss the transversality and stability, since it is 
supposed that the essential control parameter space U /G c for the continuous symmetry is identical to 
that of the discrete symmetry U /Gd in low-dimensional spacetime where the stability have been studied. 

These are class D multiunfoldings. 

5 concrete classification: reflection symmetry and axial symme- 
try 

5.1 reflection symmetry in three-dimensional spacetime 

The main purpose of the present article is to classify the stable crease set of an event horizon in an 
axially symmetric spacetime. To do so, we study Maxwell sets in reflection symmetric three-dimensional 
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spacetime where U is two dimensional and M is one-dimensional, in advance. 
A reflection isometry form a Z 2 group, that is 

tp r : M i—> M, U ^U, (46) 
G r = {id, (fi r \<Pr tpr = id} — %2, (47) 

where id is the identical map. 

5.1.1 Class A multiunfolding with reflection symmetry 

First, we consider a reflection symmetric G r -orbital multigerm of unfolding Fc T without a fixed point (a 
Class A multiunfolding). As discussed in the previous section, a G r -invariant G r -orbital multigerm of 
unfolding G r [F] is generated by the reflection isometry ip r from a usual multiunfolding F and is G r -stable 
at the minimum points if the original multiunfolding F is stable at the minimum points. 

We put the classification of multiunfolding without symmetry, that is the case of codim < 4 from 
[H] 2 , so that they are also stable in the case with two-dimensional U. A stable multiunfolding F is right 
equivalent to one of the following multiunfolding germ at their minimum points. 

• A 3 = x i + u 2 x 2 + mx, at (0, 0)6lx K 2 . 

• (j4i, j4i) = (x\ + u\, xf), at ((0i, O2), 0) G K 2 x R 2 where 0j is the origin of local coordinate (xi, 

• (Ai,Ai,Ai) = (.T 2 + Ml ,x 2 + u 2 ,x 2 ),at ((0 1; 2 , 3 ), 0) eR 3 x R 2 . 

where Ai means a stable unfolding of the stratum Ai. The unfolding A\ is omitted since its Maxwell set 
is empty set. Then, we give reflection symmetric G r -orbital multigerms of unfolding as G r (F). 
Class A multiunfolding 

A G r -invariant G r -stable multiunfolding in Class A is right SPD-equivalent to one of the following 
G r -orbital multigerms of unfolding at the minimum points. 

• G r [A 3 ] = [A 3l tp;A 3 ] = [A 3 (x,u),A 3 (ip r {x),<p r {u))]. 

. G r [{A ll A 1 )] = \{A ll A 1 ) 1 v* r {A u A 1 )] = \{A 1 ,A 1 ), {A 1 {i Pr {x),ip. r {u)),A 1 {i Pr {x), <p r {u)))}. 

. G P [(Ai,Ai,i4i)]= [(A 1) A 1 ,A 1 ), i p*(A 1 ,A 1 ,A 1 )} 
= [(A 1 ,A 1 ,A 1 ), (A 1 (<f r (x),(p r (u)),A 1 ((p r (x),(p r (u)),A 1 (<fi r (x),if r (u)))]. 

Since the function space of perturbation is same as that of without symmetry in the neighborhood of 
((pi, —,Pk), u o), these are all of the stable multiunfoldings in this class. Their Maxwell sets are given by 
the Maxwell sets in [S] and their reflectional images. They are illustrated in FigQ] 

5.1.2 Class B multiunfolding with reflection symmetry 

Next we consider Class B multiunfoldings, where uq is on a fixed point of (f r . For convenience, we 
introduce global coordinates (x,u x ,u z ) € M x U by which the reflection symmetry can be represented. 
The reflection isometry is given by (x, u x , u z ) 1— > (— x, —u Xl u z ) in our convention of the coordinates. 
We will find the list of the stable multiunfoldings in a procedure developed in the previous section. 

2 Since the dimensions of M is different, the codim is differ by one between this case and |H] 
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1. Since we treat only global minimums and M is one- and U is two-dimensional, only Ax and A3 are 
considered among the list of the strata fDef !13[l . whose codimension is less than four and the number 
of control parameter is not more than two (= dim(U)). Their stable unfoldings are A% = x 2 , and 



2. Since even after one imposes the reflection symmetry the codimension of the strata for multijet space 
Afc fl (Ai x ...) does not decrease more than [codim/2] ([x] is the Gauss's symbol), it is sufficient to 
consider the case in which codimension of A& fl (Ai x ...) is less than five (the number of control 
parameters is less than four), where the minimum value is fixed by the reflection symmetry. The 
list of the fc-tuple of the stable unfoldings is 



3. Imposing the symmetry relation (|34l) by the global coordinate, the above fc-tuples of unfoldings 
become the following G r -invariant G r -orbital multigerms of unfolding if the number of their control 
parameters ui finally is not exceeding dim(f7). (A3, Ax) is not allowed since the number of parameter 
is three even after imposing symmetry. 

4. We can concretely confirm that projection p : C°°(M x U,M) 1— > G r -C°°(M x U,M) maps neither 
the strata nor jet section into a point. Therefore the following G r -invariant G r -orbital multigerms 
of unfolding G r (F) are G r -stable at the minimum points, since their multijet section transversally 
crosses the strata Aj, fl (Ai x ...) in the jet space of G r -C°°(M x U, R). This will be proved in the 
same logic as in [Jj]. 

Class B multiunfolding 

A G r -invariant G r -stable multiunfolding in Class B is right SPD-equivalent to one of the following 
G r -orbital multigerms of unfolding at the minimum points. 

• (Ax, Ax) = (xfx + u z ,x 2 r2 ). 

To find the Maxwell set of it, we determine the minimum function fDcf ll7fl of the G r -orbital 
multigerm of unfolding, 



A3 = X 4 + U2X 2 + U\X. 



(Ax, Ax) 
(Ax, Ax, A,) 
(Ax, Ax, Ax, Ax) 

(40 

(A 3 ,Ax) 



(x\ +ui,ac|) 

I 2 2 1 2\ 

(X x + Ux,X 2 + U 2 , X 3 ) 

I 2 2 1 2 1 2\ 

(x x + ux, x 2 + u 2 , x 3 + u 3 , x 4 ) 

(x\ + u 2 x\ + UxXx) 

(x\ + u 2 xf + UxXx,x\ + U3). 



(48) 
(49) 
(50) 
(51) 
(52) 




u z < 
u z >0 



(53) 



Then the Maxwell set that is the singular point of the minimum function is 



-^Maxwell (Ax, Ax) = {(u x ,u z )\u z 



0}. 



(54) 



• G r (Ax) = (Ax, Lp* r Ax) = ((x r - x ) 2 + u x , (x r + x ) 2 



u x ) ~ (x\ + u x , x\ — u x ), 



where x\ and x 2 are local coordinates. 
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The minimum function is 



Tam{Ax,<p r Ax) = < (55) 
u x u x > 



The Maxwell set that is the singular point of the minimum function is 

-B M axwcii(^i, ¥>tAi) = {(u x ,u z )\u x = 0}. (56) 



• (G r (Ai),Ai) = ((Ai,(p*Ai), Ai) = (x 2 +u x + u z ,x 2 - u x + u z ,x 2 3 ). 



The minimum function is 

Ta\n{{Ai,ip* r Ai),Ai) = min(-\u x \ +u z ,Q) = \ ^ j Ma: j ^ Uz (57) 

[0 \u x \ < u z 

The Maxwell set that is the singular point of the minimum function is 

fiMaxweu((Al,^*Al), Al) = {{u x ,U z )\u x = 0, U z < 0} 

U {{u x ,u z )\\u x \ = u z }. (58) 
. G r {{A t , A x )) = {G r {A\), G r {Ai)) = {{At, <p*Ai), {A lt 

= ((^i + U x , X2 — "1), (^3 + u x + u z , x^ u x + w z )), 
where xx,X2,X3, x± are local coordinates. 

The minimum function is 

I Z r, ■ ( 59 ) 

-\u x \ \u z > 

The Maxwell set that is the singular point of minimum function is 

-BMaxwcll((Al,V ; '*^l)7(^l:^*^l)) = {{u Xl U z )\u x = 0} 

U{{u x ,u z )\u z =0}. (60) 

• G r (A 3 ) = A 3 = x^ + u z x 2 + u x x r . 

The Maxwell set is same as that of A 3 without symmetry. 

-B M axwcii(^3) = {{u x ,u z )\{u z <0,u x = 0)}. (61) 

It should be noted that some multiunfoldings decreases the number of control parameter after imposing 
the reflection symmetry. Especially it is note worthy that {{Ai, p*Ai), {Ai, (p*A\)) is not stable without 
the reflection symmetry but stable with the reflection symmetry. Then the entirely new Maxwell set have 
been added to the classification. These Maxwell sets are illustrated in Fig|2] 

5.2 axial symmetry in four-dimensional spacetime 

Here we consider axial symmetry in a four dimensional spacetime. In this case, its isometry group G a is 
continuous group and isomorphic to U{1), 

d 



(fg : M >—¥ M, U 1— > U, ipg = exp 
G a = {<p e \0<6<2Tr,<p 2w = id}~U{l), (62) 
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5.2.1 Class C multiunfolding with axial symmetry 

For the reflection symmetry, the function space was restricted by only finite codimensions and there 
does not appear an essentially new unfolding without a fixed point. As the axial symmetry, however, is 
continuous, our standard investigation will fail by two reasons. Since the function space (strictly speaking, 
its jet space) is restricted into infinite codimensional subspace, a new function which never provides stable 
unfolding owing to its infinite codimension can provide a stable unfolding, e.g. (x 2 + y 2 ) 2 see To find 
all of these new functions it is required to survey the whole right SPD equivalence class of the restricted 
function space systematically and it is embarrassing. 

Another is that infinite number of stationary points can be generated by the axial isometry as ip* a F '. 
Then we cannot use semi-local discussion where usual function germ defined. Nevertheless, the essential 
control space of this axial symmetry in four dimensional spacetime U/G a is identical to that of reflection 
symmetry in three-dimensional spacetime U' /G r . Omitting the angular direction of the state space M 
and control space U , we treat the remaining as a state space and control space U/G a . We can study it 
in semi-local formalism. 

To determine a correspondence between U r = (U' ~ M. 2 )/G r and U a — (U ~ M. 3 )/G a , it is convenient 
to use coordinate systems (u' x , u' z ), [u r , ug,u z ). By identifications u' x = u r and u' z = u z , a correspondence 
K 2 /G r I— > R 3 /G a is given, since the conditions < u' x < oo, — oo < u' z < oo are consistent to < u r < 
oo, — oo < u z < oo. This correspondence simply gives the correspondence between Class A multiunfolding 
and the Class C multiunfolding. We also use local coordinate (xi,Ui) to show local diffeomorphism 
equivalence, xi and u, are local coordinate in r-z or u r -u z section of M and U, respectively. 

To determine a Maxwell set we should give a unfolding in U not in U/G a . From the splitting 
lemma ITS], we understand that it is sufficient to add a Morse function by a dummy variable ug of angular 
direction, which does not concern the structure of the minimum point, as SO 2 = (0 — ug) 2 where — ug 
is considered to be small 3 . 

Class C multiunfolding 

A G a -invariant G a -stable multiunfolding in Class C is right SPD-equivalent to one of the following 
G a -orbital multigerm of the unfolding at the minimum points. 

• G a [(A lt Ai)] = (x 2 +y 2 1 +u 1 + {6 l - ugf,x 2 +y 2 + (9 2 - ug) 2 ), 
< 0i < 2vr, < 9 2 < 2tt, < ug < 2vr. 

For G a [(Ai, Ai)], the minimum function is 

mm(G a [(A 1 ,A 1 ))) = T 1 Ul <° . (63) 
10 ui > 

Then its Maxwell set is given by 

Bua,xwen(Ga[(Al,Ai)]) = {(ui , u 2 , ug) \m =0}. (64) 

• G a [{Ai, A u Ai)} = {xj +y\+ui + {6 X - u g ) 2 , x 2 +y 2 +u 2 + {6 2 - u e ) 2 ,x 2 +y 2 + (6» 3 - u e ) 2 ), 
< 0i < 2?r, < 6 2 < 2tt, < 6 3 < 2?r, < ug < 2n. 

3 One should not translate this as the generator of an event horizon does not change its angular direction. To write this 
form, we may change the origin or scale of the angular coordinate 8 
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For G a [(Ai, Ai, Ai)], the minimum function is 

{U\ U\ < 1*2, Ui < 
u 2 u 2 < u x ,u 2 < Q ■ (65) 
ui > 0, u 2 > 

This is singular at «i — u 2 < 0, ui = 0, u 2 > and u 2 = 0, i*i > 0. Then its Maxwell set is given 

by 

BMaxwen(G a [(Ai, Ai, Ai)}) = {(ui,u 2 ,ue)\ui = u 2 < 0} 

U {(u-L,u 2 ,ug)\ui — 0,u 2 > 0} U {(ui,U2,ug)\u 2 — 0, u x > 0}. (66) 

• G a [A 3 ] = x\ + u 2 x\ + mxi + (9i - u e ) 2 , 
< 0i < 2tt, < u g < 2tt. 

For G a [A3], the Maxwell set is same as that of A3 on the section ug = 0. Then its Maxwell set is 
given by 

SMaxwcii(G a [74 3 ]) = {(ui,u 2 ,ug)\u 2 > 0,iti = 0} (67) 
These Maxwell sets are illustrated in Fig|3] 
5.2.2 Class D multiunfolding with axial symmetry 

Similarly to the Class C multiunfolding, the correspondence between U r — M 2 /G r and U a = M 3 /G a also 
gives correspondence of multiunfolding between Class B and Class D. In addition that, it is necessary 
to show that the minimum is not isolated but make an orbit of the isometry. This is realized by global 
coordinates of the symmetry so that the Morse function S9 2 , which is added to revive a dummy variable 
ug, is extended to a global function of 89. It is achieved by cos SO ~ 1 — S9 2 /2 + .... 
Class D multiunfolding 

A G a -invariant G a -stable multiunfolding in Class D is right SPD-equivalent to one of the following 
G a -orbital multigerm of unfolding at the minimum points. 

• {Ax,Ai) = (r 2 +u z ,r 2 ). 

For (Ax,Ai), the minimum function is 

min(A u A 1 ) = h Uz <°. (68) 
10 u z > 

The Maxwell set that is the singular point of minimum function is 

BMaxwellifeAi) = {(u r , u z , ug)\u z = 0}. (69) 

The toroidal event horizon reported in ref. |5j is caused by the crease set which is embedding of this 
Maxwell set. 

• G a (Ai) = z 2 + u r cos(/9 — ug), 
< 6><2tt, < ug < 2ir. 
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For G a (Ai), the minimum function is 



min(G B (i4i))=u r . (70) 

This is singular at u r = 0, 

SMaxwcii(G Q (Ai)) = {(u r ,u z ,u 9 )\u r = 0}. (71) 

This Maxwell set is reported in the collision of the two black holes 

• (G a (Ai),A!) = {z 2 + u r cos(9 -u e ),rl+u z ), 
< 6><2tt, < lie < 2?r 

The minimum function of (G a (Ax), A\) is 

min(G a (A 1 ), A x ) = min(u r ,u z ) = \ " 2 Ur Uz . (72) 

I u r u r < u z 

This is singular at u r = 0, u z > and u z = u r ,u z < 0; 

Buaxwell(G a (Ai)) = {{u r ,U z , Ue)\u r = 0, U z > 0} 

U {(u,.,u z ,u e )|u z < 0,u z = u r }. (73) 

• (G a (Ai), G a (A 1 )) = (z\ + u r cos(6»i - u e ) +u z ,z 2 + u r cos(6* 2 - u e )), 
< 6\ < 2?r, < 9 2 < 2tt, < u e < 2tt. 

The minimum function of (G a (Ai), G a (Ai)) is 

min((G„( J 4i) ) G (Ai)))=ndn(u r + u, s « r )= ^ r * . (74) 



it r + u z < u z 



This is singular at u r — and u z = 0; 



S M axwell((G a (Ai),G a (Ai))) = {(u r ,u 2 ,w e )|u r = 0} 

U {{u r ,u z ,ug)\u z = 0}. (75) 

• G a (A 3 ) = A 3 = r 4 + u z r 2 + u.r, 

where r and u is the position vector to (r, 9) and {u r , ug) surface of M and U, respectively. 

For G a (A 3 ), the Maxwell set is same as that of A 3 

-BMaxweii (G a (A 3 )) = {(u r ,u z ,ug)\u r — 0,u z < 0}. (76) 

Their Maxwell sets are illustrated in FigQ] 

Here it should be noted that the Maxwell set stably possesses one-dimensional component G a {A\), 
(G a (Ai), Ai), (G a (Ai), G a (Ai)) or G a (A 3 ) 1 while only two-dimensional structure is stable without symmetry [5] 0. 

From the above classifications, we are able to see the crease set. As stated in the crease set is 

an acausal embedding of the Maxwell set into the spacetime so that the crease set is homeomorphic to a 
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ball and approaches to null at its boundary. We give a figure of an example of axially symmetric Maxwell 
set in Fig|SJ 

In the example in Fig|Sl by an appropriate timcslicing, it is observed that toroidal black hole and 
spherical black hole coalesce. At the lower spatial hypersurface, the section of the crease set is a ring and 
a point. They would be a crease of a toroidal event horizon and a cusp on a spherical event horizon. As 
time goes by they approach each other and the ring shrinks, and then they unite and disappear. This 
implies the coalescence of toroidal black hole and a spherical black hole. 

To find a example in Fig03 we have combined the structures of the Maxwell set so that a total 
Maxwell set is homeomorphic to a ball, since it is required by the assumption that the event horizon 
finally becomes a single sphere. Nevertheless it might be important to consider topologically non-trivial 
final state of the event horizon, which becomes meaningful in studies of higher-dimensional black hole. 
For example, for a toroidal final event horizon, we compose the Maxwell set so that it is homeomorphic 
to a ring. 

6 conclusion and discussions 

We have investigated the Maxwell set under a discrete symmetry and a continuous symmetry. In the 
context of singularity theory, the stable Maxwell set is classified. Then we concretely showed the clas- 
sification of the stable Maxwell set with reflection symmetry in three dimensional spacetime and axial 
symmetry in four-dimensional spacetime. Remarkable fact is that by the axial symmetry one-dimensional 
segments of the crease set is possible to become stable while it is not stable in four-dimensional spacetime 
without symmetry. 

It might be worth to comment that the present result also is relevant in the case of almost symmetric 
where the symmetry is dynamically stable, if a very fine structure is neglected. Some cases will be almost 
axially symmetric in realistic gravitational collapse. Nevertheless it is not mathematically clear how can 
we neglect the fine structure. Under the present situation, our result only suggests the structure of the 
crease set in almost symmetric spacetime. 

By the way, one may be doubtful that peculiar crease sets for example like in Fig|S]can be formed. 
However, it will be realistic in the collision of collapsing stars. In an oblate spheroidal event horizon, it is 
probable that the crease set is naturally two dimensional[5][H]- If sufficiently oblate spheroidal collapsing 
star and spherical collapsing star are sufficiently away in direction of their symmetry axis and if they 
coalesce, the crease set like FigE^will be formed. Furthermore, when several black holes coalesce, their 
crease set will become fairly complicated. 

Since these collapses are axially symmetric, it will be possible to simulate them by numerical cal- 
culation of the gravitational collapse. Not only the fairly simple case in refs.00 (these cases are also 
reflection symmetric about x-y plane) but also some considerably complex cases should be studied. With 
that the various topological black hole will be observed simultaneously. That will be thought-provoking 
for the gravitational collapse and phenomena concerning the event horizon, e.g. quasi-normal mode of 
gravitational radiation. 
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(G r (Aj), G r (Aj)) 



Figure 2: The Maxwell sets in two-dimensional reflection symmetric spacetime on its fixed points. They 
are the Maxwell sets of Class B multiunfoldings. 





i 



24 




25 



(A,, a,) m 



G a (A,) 



G a (A 3 ) 



(G a (.AMi) 




(GJA,), G a (A,)) 




Figure 4: The Maxwell sets in three-dimensional axially symmetric spacetime on its fixed points. They 
are the Maxwell sets of Class D multiunfoldings. 
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Figure 5: An example of the axially symmetric Maxwell set is shown. That is composed of G a {A^), 
(G a (Ai), Ai) and G a [^3], which are connected by G a {A\) and G a L4i, A{\. Next figure is the crease set 
which is given by embedding the Maxwell set into a spacetime as a spatial set. The bottom figure indicate 
the change of the topology of an event horizon with that crease set. 
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